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We  developed  mathematical  models  for  problems  involving  the  fabrication  of  ceramic- 
matrix  composites,  the  coating  of  ceramic  fibers,  and  the  oxidation  of  ceramic-matrix 
composites.  To  solve  these  models,  we  developed  an  adaptive  and  parallel  computational 
framework  for  the  finite  element  and  finite  volume  solution  of  partial  differential 
equations.  The  parallel  adaptive  framework  handles  mesh  generation,  refinement,  and 
coarsening;  dynamic  load  balancing;  distributed  database  management;  interprocessor 
communications;  data  migration;  time  integration;  and  algebraic  solutions.  The 
models  and  solutions  revealed  improved  processing  techniques  for  a  reactive  vapor 
infiltration  application  and  potentially  superior  strategies  for  coating  fibers. 
Oxidation  studies  are  investigating  the  "self  healing"  of  cracks. 
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1.  Objectives 

Goal:  The  twofold  objectives  of  this  project  were  (i)  modeling  and  analysis  of  prob¬ 
lems  in  ceramic  material  fabrication  and  (ii)  the  development  of  adaptive  and  parallel  tech¬ 
niques  to  automate  the  solution  of  these  problems. 

Problems  involved  (i)  fabrication  of  ceramic-matrix  composites  by  vapor  infiltration 
and  deposition  processes,  (ii)  the  coating  of  ceramic  fibers  by  vapor  deposition,  and  (iii) 
estimates  of  the  longevity  of  ceramic  composites  when  subject  to  oxidation  at  elevated 
temperatures.  Although  specific  processing  problems  were  addressed,  the  models  and 
parallel  adaptive  software  have  a  much  wider  range  of  applicability.  The  adaptive  tech¬ 
niques  combine  mesh  refinement  and  coarsening  (h-refinement),  method-order  variation 
(p-refinement),  and  mesh  motion  (r-refinement)  in  both  space  and  time.  The  software  exe¬ 
cutes  on  seritd  and  parallel  computers  and  operates  under  the  Message  Passing  Interface 
(MPI).  It  is,  thus,  portable  and  efficient. 

Highlights  of  our  findings  are  presented  in  Section  2  and  more  detailed  descriptions 
follow  in  Section  3.  Supported  personnel  are  listed  in  Section  4,  Publications  appear  in 
Section  5,  and  presentations,  interactions,  and  transitions  appear  in  Section  6. 

2.  Significant  Accomplishments  and  Findings 
2.1.  Parallel  Adaptive  Computation 

Adaptive  computation  provides  the  most  effective  way  of  generating  reliable, 
efficient,  and  robust  solutions  of  large-scale  problems  involving  partial  differential  equa¬ 
tions.  Beginning  with  a  trial  solution  generated  on  a  coarse  mesh  with  a  low-order 
method,  an  adaptive  computation  appraises  the  accuracy  of  a  solution,  enriches  the  solu¬ 
tion  where  needed  by  combinations  of  h-,  p-,  and  r-refinement,  and  recursively  generates 
solutions  until  prescribed  accuracy  requirements  have  been  satisfied.  Error  estimates  con¬ 
trol  adaptive  enrichment  and  provide  a  requisite  measure  of  solution  accuracy.  Without 
such  error  estimates,  computations  must  either  be  repeated  with,  e.g.,  finer  meshes  or  solu¬ 
tion  accuracy  must  be  inferred  using  independent  knowledge  of  the  anticipated  solution. 
The  former  approach  is  inefficient  and  impractical  with  truly  demanding  problems  while 
the  latter  approach  provides  qualitative  information  at  best  and  is  becoming  less-and-less 
possible  with  the  complexities  of  modem  engineering  systems.  Error  estimates  obtained 
by  order  embedding  (p-refinement)  are  asymptotically  correct  and  have  expenses  ranging 
from  one-third  to  one-half  of  the  solution  cost.  Combinations  of  h-  and  p-refinement  are 
remarkably  efficient  and  may  have  exponential  convergence  rates. 

Parallel  computation  is  essential  for  demanding  multi-dimensional,  nonlinear  prob¬ 
lems.  However,  developing  parallel  adaptive  software  is  quite  a  challenge.  Adaptive  pro¬ 
cedures  use  nonuniform  stmctures  and  complex  heuristics  to  achieve  efficiencies.  These 
are  difficult  to  parallelize  using  standard  constmcts.  Load  balancing  must  be  dynamic, 
since  adaptive  enrichment  will  upset  a  balanced  computation 

We  have  developed  a  parallel  adaptive  computational  framework  for  solving  three- 
dimensional  transient  and  steady  partial  different!^  equations  by  finite  element  and  finite 
volume  schemes.  The  software  has  overcome  many  of  the  cited  difficulties  and  contains 
several  novel  and  innovative  features  including: 
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•  Automatic  three-dimensional  grid  generation  by  recursive  octree  decomposition  with 
an  interfaces  to  several  CAD  systems  [11].* 

•  Adaptive  h-,  p-,  and  r-refinement  procedures  [11-13].  (Adaptive  p-refinement  is  not 
yet  operation^  in  three  dimensions.) 

•  A  mesh  database  for  the  efficient  storage  of  information  associated  with  high-order 
finite  element  computation  [11]. 

•=  Dynamic  load  balancing  by  iterative  tree  [3-4,  7-8,  10-13],  inertial  bisection  [11],  and 
octree  bisection  [7]  procedures. 

•  A  parallel  mesh  database  to  handle  the  migration  of  data  and  communications 
between  processors  [4,  8,  11]. 

•  Adaptive  h-  and  p-refinement  procedures  for  temporal  integration  that  include 
backward-difference  and  singly  implicit  Runge-Kutta  methods. 

•  Computational  linear  algebra  routines  that  interface  with  the  mesh  and  parallel  mesh 
databases. 

•  Visualization  software  using  Data  Explorer. 

This  is  the  first  adaptive  system  containing  these  capabilities.  The  framework  uses 
the  object-oriented  methodology  to  ease  development  and  maintenance  and,  as  noted,  the 
MPI  for  portability.  Scalability  has  been  verified  by  extensive  computation  on  fluid  flow 
[12]  and  crystal  growth  problems.  Although  this  grant  has  ended,  work  on  this  effort  is 
continuing  as  part  of  an  AFOSR  MURI  grant.  Additional  development  will  enhance  the 
time  integration  and  linear  algebraic  procedures,  improve  the  migration  strategies,  and 
optimize  the  hp-  and  hpr-refinement  techniques.  The  final  framework  will  simplify  the 
development  of  new  finite  element  applications.  Many  will  only  require  software  for  the 
generation  of  element  stiffness  matrices  and  preconditioners. 

2.2.  Composite  Materials  Fabrication 

Working  with  colleagues  at  Rensselaer’s  Materials  Science  and  Engineering  Depart¬ 
ment,  the  Air  Force’s  Wright  Laboratories,  Centric  Engineering,  and  Dupont  Lanxide,  we 
modeled  processes  for  manufacturing  ceramic  matrix  composites,  coating  ceramic  fibers, 
and  oxidizing  ceramic  composites  at  elevated  temperatures. 

Ceramic  composites  offer  several  advantages  relative  to  homogeneous  materials  when 
used  in  high-temperature  applications  such  as  heat  exchangers,  heat  engine  components, 
and  aerodynamic  surfaces  exposed  to  hypersonic  flows.  Their  fabrication  by  a  variety  of 
techniques  (e.g.,  plasma  spraying,  sintering,  or  chemical  vapor  infiltration)  remains 
difficult  Processing  times  are  long  and  production  costs  are  high.  Mathematical  model¬ 
ing  and  computational  experimentation  are  identifying  improved  or  optimal  materials  and 
fabrication  techniques  for  military  and  civil  applications. 

As  noted,  ceramic  composities  operate  in  hostile  environments  and  are  subject  to 
chemical,  thermal,  and  mechanical  attack.  Bonds  between  fibers  and  the  matrix  can  fail. 
Oxidizing  reactions  can  weaken  or  destroy  both  fiber  and  matrix.  Fiber  coatings  can 


*  Cited  references  appear  in  Section  5. 
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provide  better  adhesion  with  matrix  materials  as  well  as  protect  fibers  from  chemical 
attack.  Cracks  in  the  matrix  can  destroy  protective  fiber  and  matrix  coatings  and  expose 
the  interior  of  the  composite  to  the  caustic  environment.  Materials  may  oxidize  during 
prolonged  exposure  and  this  may  seriously  degrade  the  mechanical  properties  of  the  com¬ 
posite.  Alternatively,  volume  expansion  accompanying  reaction  may  close  (or  “heal”)  the 
crack  before  substantial  damage  can  occur.  Again,  mathematical  modeling  can  illuminate 
these  phenomena  and  suggest  materials  and  strategies  that  limit  adverse  effects. 

We  have  investigated  these  issues  and  our  key  findings  and  results  follow. 

•  We  developed  a  model  for  fabrication  of  ceramic  composites  by  reactive  vapor 
infiltration  (RVI)  that  couples  reaction  with  mechanical  deformation.  The  model  and 
adaptive  software  have  predicted  pore  closing  between  powder  grains,  deformation 
caused  by  reaction,  and  residual  stresses  [6,  9,  15]. 

•  Using  the  RVI  model,  we  showed  that  fabrication  times  and  excessive  deformation  of 
molybdenum-disilicide-matrix  composites  may  be  reduced  by  initiating  the  process 
with  a  mixture  of  50%  powdered  molybdenum  and  50%  molybdenum  disilicide. 
This  strategy  was  bom  out  experimentally  and  the  resulting  composites  had  excellent 
mechanical  properties  [9,  15]. 

•  We  developed  a  model  for  the  coating  of  ceramic  fibers  by  chemical  vapor  deposition 
(CVD).  The  model  has  capabilities  to  analyze  transient  and  steady  flows  containing 
multiple  species;  chemical  surface  reactions;  heating  in  hot-  and  cold-walled  reactors 
by  conduction,  convection,  and  radiation;  and  geometric  variation  due  to  fiber  motion 
and  deposition  [17]. 

•  The  fiber  coating  model  has  been  used  to  study  the  coating  of  sapphire  fibers  with 
beta-alumina  in  hot-walled  reactors.  We  examined  critical  parameters  (flow  rate, 
fiber  speed,  precursor  injection  rates,  reactor  temperature  and  pressure)  and  suggest 
strategies  that  minimize  losses.  Since  the  walls  of  the  reactor  are  heated  with  the 
fiber,  they  too  are  coated.  We  developed  parameter  choices  that  minimize  wall  coat¬ 
ing  and  process  time  while  producing  a  uniform  fiber  coating  [17]. 

•  Researchers  at  the  Air  Force’s  Wright  Laboratories  are  using  our  fiber  coating  model 
and  software  to  design  a  control  system  for  their  experimental  facility.  They  will  use 
the  model’s  predictions  to  create  a  neural-network  model  of  the  parameter  space  that 
will  provide  an  active  control  strategy  during  processing. 

•  We  developed  a  model  for  the  oxidation  of  ceramic  composites  at  elevated  tempera¬ 
tures.  The  model  addresses  the  flow  of  oxidants  in  a  crack  or  pore  and  their  reaction 
with  and  diffusion  through  surfaces  to  oxidize  the  matrix  material.  It  has  capabilities 
to  handle  an  arbitrary  number  of  chemical  species,  surface  and  bulk  reactions,  and 
fluid  and  solid-state  diffusion.  It  is  coupled  with  a  mechanical  deformation  model  to 
predict  stress  concentration  and  further  cracking  or  the  healing  of  cracks  by  expan¬ 
sion  of  the  matrix  [18]. 

•  The  model  is  being  used  to  study  the  oxidation  of  silicon-carbide  matrix  composites 
by  oxygen  and  water  vapor.  This  work  is  in  collaboration  with  Materials  Scientists 
at  Rensselaer  and  at  Dupont  Lanxide. 

All  models  were  discretized  and  solved  by  finite  element  methods  using  the  adaptive 
and  parallel  framework.  The  teaming  of  this  theoretical  analysis  with  experimental 
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investigation  enabled  knowledge  to  obtained  more  rapidly  than  each  would  have  alone. 
The  theoretical  studies  identified  parameter  values  that  were  likely  to  lead  to  improve¬ 
ments.  These  computational  investigations  took  far  less  time  than  experimentation,  which 
often  requires  several  days.  Experimentation,  on  the  other  hand,  provided  the  necessary 
data  and  constitutive  parameters  to  the  mathematical  model. 

3.  Models  and  Methods 

3.1.  Adaptive  Computation 

The  adaptive  and  parallel  computational  framework  for  partial  differential  systems 
has  the  generic  form 

a,m(x,r,q)  -i-/(x,t,q,Vq)  =  V-d(x,r,q,Vq) 

where  x  and  t  denote  space  and  time,  q(x,r )  denotes  the  solution  vector,  8,  denotes  partial 
differentiation  vdth  respect  to  t ,  and  V  denotes  the  gradient  and  divergence  operators. 
The  mass  rnatrix”  8qm  and  the  “diffusion  matrix”  8vqd  may  be  singular  so  the  system 
may  be  of  mixed  type.  The  above  equations  are  discretized  by  the  finite  element  method 
in  space  and  by  either  singly-implicit  Runge-Kutta  (SIRK)  methods  or  backward 
difference  formulas  (BDFs)  in  time.  These  combinations  have  been  successfully  used  to 
solve  problems  involving  fluid  flow  [3,  8,  10,  11,  13]  and  crystal  growth  in  addition  to  the 
materials  piwessing  problems  described  herein.  High-order  finite  element  bases  are  avail¬ 
able  using  hierarchical  bases  where  polynomial  representations  of  the  solution  are  regarded 
as  corrections  to  lower  order  polynomials. 

Special  stabilization  terms  are  available  to  handle  systems  where  reaction  and/or  con¬ 
vection  dominate  diffusion.  This  is  the  case  with  the  materials  fabrication  problems  under 
investigation.  When  not  done,  solutions  of  these  systems  exhibit  spurious  oscillations. 
The  stabilization  model  used  for  this  study  was  developed  in  collaboration  with  research¬ 
ers  at  Centric  Engineering.  As  noted,  SIRK  and  BDF  techniques  are  used  for  time 
integration.  These  may  be  done  as  a  method  of  lines  (MOL)  formulation,  where  the 
method  and  time  step  are  independent  of  the  spatial  location,  or  as  a  local  refinement 
meAod  (LRM),  where  different  methods  and  time  steps  are  used  in  different  spatial 
regions.  The  latter  alternative  is  being  implemented  in  three  dimensions,  but  its  perfor¬ 
mance  in  two  dimensions  has  shown  it  to  be  a  very  effective  method.  SIRKs  have 
remarkable  stability  properties  and  efficiencies  rivaling  those  of  BDFs.  They  shine  when 
used  with  LRMs  and  in  evolving  geometry  problems. 

Estimates  of  spatial  discretization  involve  the  comparison  of  solutions  obtained  with 
different  order  finite  element  approximations.  Localization,  made  possible  through  super¬ 
convergence,  simplifies  the  computation  and  enables  estimation  at  the  element  level  [5, 
1^3*  Accuracy  appraisal  is  used  to  drive  the  adaptive  procedure,  but  is  also  important  to 
provide  assurance  that  the  solution  is  being  calculated  correctly. 

Algebraic  solutions  are  obtained  using  direct  methods  (sparse  Gaussian  elimination) 
for  smaller  problems  and  generalized  residual  methods  with  preconditioning  for  larger 
problems.  Preconditionings  are  difficult  to  generalize  and  currently  involve  multilevel  and 
relaxation  techniques. 
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3.2.  Parallel  Computation 

To  operate  effectively  in  a  parallel  solution  environment  each  of  the  components  of 
an  automated  adaptive  analysis  system  must  operate  in  parallel  and  maintain  computational 
effectiveness  as  the  problem  scales.  Key  to  maintaining  parallel  performance  is  the  inclu¬ 
sion  of  procedures  that  directly  account  for  the  adaptive  nature  of  the  numerical  discretiM- 
tion,  and  that  maintain  the  appropriate  balance  of  computational  work  load  over  the  pro¬ 
cessors. 

Current  load  balancing  methods  involve  iterative  neighboriiood  relaxation  [3-4,  7-8, 
10-13],  inertial  bisection  [11],  and  octree  bisection  [7].  All  execute  in  dynamic  environ¬ 
ments  that  restore  a  balanced  computation  that  was  interrupted  by  adaptivity.  The  former 
constructs  and  balances  load-request  trees  where  each  processor  requests  load  fix)m  its 
most  heavily  loaded  neighbor.  The  latter  two  procedures  recursively  bisect  the  mesh  in 
directions  perpendicular  to  its  principal  axis  of  inertia  and  coordinate  axes,  respectively. 
Mesh  and  solution  data  are  migrated  between  processors  by  parallel  sort  and 
communication-reduction  schemes  [4,  1 1]. 

The  framework  has  capabilities  to  support  parallel  hp-  and  hpr-refinement  methodolo¬ 
gies  for  multi-dimensional  problems.  Future  versions  will  use  object-oriented  methodol¬ 
ogy  to  provide  commonality  and  reuse  by  abstracting  many  operations  and  algorithms  to 
allow  them  to  be  combined  in  various  ways  to  produce  an  application.  Examples  of  com¬ 
mon  items  are  time  integration  procedures,  linear  and  nonlinear  algebraic  solution  tech¬ 
niques,  spatial  integration  routines,  mapping  procedures  from  physical  to  the  parametric 
space  of  an  element,  and  interpolations  of  field  variables  over  the  mesh  entities  (volumes, 
faces,  edges,  and  vertices). 

When  the  geometric  domain  evolves  with  the  solution,  its  representation  must  be 
properly  updated  to  account  for  appropriate  smoothness  requirements  and  surface  motion. 
Basic  tools  to  update  smooth  geometric  representations  as  dictated  by  the  motion  of  the 
discrete  mesh  representation  are  available;  however,  they  are  limited  in  their  ability  to 
effectively  support  irregular  meshes  and  computational  procedures.  Development  will  be 
necessary  to  update  geometric  model  surface  representations  to  support  remeshing  or  other 
operations  needing  geometric  information  such  as  curvature. 

3.3.  Reactive  Vapor  Infiltration 

With  Rensselaer  Materials  Scientists  William  Hillig  and  John  Hudson,  we  studied 
RVI  techniques  for  manufacturing  ceramic  composites  where  silicon  carbide  (S/C)  or 
alumina  (A/2O3)  fibers  are  mixed  with  molybdenum  (Mo)  powder  or  a  mixture  of  Mo  and 
molybdenum  disilicide  (MoSi^  and  pressed  into  a  porous  preform.  The  preform  is 
exposed  to  a  silicon  tetra-chloride  (S/C/4)  hydrogen  flow  where  molecular- 
surface  reactions  liberate  Si  which,  when  absorbed  into  the  preform,  reacts  with  the 
powder  to  form  the  MoSii  matrix.  We  developed  a  model  of  this  process  that  accounted 
for  the  conservation  of  mass  of  the  various  chemical  species,  the  diffusion  of  Si  into  the 
porous  media,  the  formation  of  an  intermediate  silicide  Mo^Si^,  and  the  deformation  of 
the  viscous  media  that  results  from  pores  filling  between  powder  grains  [6,  9,  15]. 

When  an  initial  Mo  powder  was  compressed  to  a  porosity  of  45%,  the  siliciding 
reactions  produced  a  158%  volume  increase  that  filled  pores,  but  was  large  enough  to 
swell  the  material  and  cause  cracking.  We  discovered  that  pores  could  be  filled  with 
essentially  no  additional  material  expansion  and  that  production  times  could  be  cut  in  half 
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by  initiating  the  process  with  a  powdered  mixture  of  50%  Mo  and  50%  MoSi2  [9]. 

Some  of  our  findings  [9]  are  illustrated  in  Figure  1.  At  the  upper  left,  we  present  a 
micrograph  of  a  cut  section  through  a  partially  silicided  pellet.  The  top  and  right  portions 
of  this  micrograph  contain  the  desired  MoSii  matrix;  the  left  and  bottom  portions  contains 
unreacted  Mo;  and  the  band  running  from  the  left  edge  to  the  bottom  contains  Mo^Si^. 
The  hardness  indentations  show  the  increased  density  in  the  MoSi2  portion  of  the  speci¬ 
men.  The  computed  solution  at  a  comparable  time  is  shown  at  the  upper  right  of  Figure 
1.  The  color  red  indicates  M0S12,  yellow  indicates  Mo^Si^,  and  blue  indicates  unreacted 
Mo  powder.  The  computational  mesh  is  finer  near  the  reaction  front  (Mo^Si^)  and  moves 
inward  (towards  the  bottom  left  comer)  as  the  reaction  progresses.  Material  distortion  is 
apparent  in  mesh  lines  that  were  originally  parallel  to  the  coordinate  axes.  Computed 
results  were  excellent  and  volumes  and  thicknesses  of  the  M0S12  and  M05S/3  layers  as 
functions  of  time  differed  by  less  than  10%  [6,  9]  from  experimental  observations. 

In  the  lower  portion  of  Figure  1,  we  display  a  pore  between  two  grains  of  Mo 
powder  at  two  times.  The  pore  is  closing  as  the  reaction  proceeds  (left  to  right)  and 
closes  completely  after  approximately  two  seconds  on  the  sUiciding  reaction.  As  pores 
close,  siliciding  proceeds  by  solid-state  diffusion.  This  is  inefficient  and  a  method  to  keep 
the  pores  open  is  needed.  This  might  be  possible  by  using  microwave  heating. 

3.4.  Ceramic  Fiber  Coating 

Researchers  in  Steven  Le  Clair’s  manufacturing  group  at  Wright  Laboratory  have 
been  studying  the  coating  of  continuous  sapphire  fibers  by  beta-alumina  using  CVD  in 
hot-walled  reactors.  We  formulated  a  model  for  their  process  using  the  Navier-Stokes  and 
energy  equations  to  model  the  flow  of  carrier  gas  through  the  reactor.  Precursor  species 
that  react  on  heated  surfaces  to  form  the  coating  are  passively  convected  by  the  carrier 
gas.  The  surface  reaction  model  includes  capabilities  for  molecular  surface  impingements 
with  and  without  successful  reactions.  This  system  is  coupled  with  an  energy  equation 
and  thickness  variation  model  for  the  fiber.  The  fiber  may  be  heated  by  conduction,  con¬ 
vection,  and/or  radiation.  As  noted,  the  Navier-Stokes  equations  have  been  stabilized  by 
the  addition  of  a  least-squares  local  diffusion  model  that  prevents  spurious  oscillations 
without  introducing  excessive  diffusion.  This  latter  aspect  of  our  investigation  was  done 
in  collaboration  with  Centric  Engineering. 

While  experimental  results  at  Wright  Laboratories  are  not  yet  available,  the  model 
and  its  computational  results  appear  to  be  qualitatively  correct.  Parameter  values,  within 
experimental  range,  have  identified  strategies  that  reduce  wall  coating  while  producing 
filers  having  uniform  coatings.  Concentrating  the  precursor  near  the  fiber  and  away  from 
the  hot  wall  provides  the  most  substantial  in^rovement  in  the  efficiency  of  hot-walled 
CVD  reactors.  Higher  reaction  probabilities  yield  higher  efficiencies,  especially  when  the 
precursor  enters  near  the  fiber.  Longer  dwell  times  (the  time  that  the  preciusor  remains  in 
the  reactor),  controlled  by  either  the  reactor  length  or  flow  rate,  improve  production  rates 
but  reduce  efficiency. 

Results  in  Figure  2  show  the  concentration  of  the  precursor  as  a  function  of  Froude 
number 

Fr  =  iPligL) 


-  8  - 


where  U  is  sl  reference  velocity,  g  is  the  acceleration  of  gravity,  and  L  is  the  reactor 
length.  In  each  drawing,  the  fiber  is  located  on  the  left  and  the  wall  is  located  on  the  right 
of  the  cylindrical  reactor.  The  illustrated  reactor  has  a  cold  wall  with  the  fiber  heated  by 
resistive  heating.  The  Reynolds  number  (Re),  Prandtl  number  iPr),  and  flow  rate  (5c)  are 
indicated  at  the  bottom  of  the  figure.  Other  parameters  are  recorded  in  Adjerid  et  al.  [17]. 
An  axisymmetric  steady  solution  was  generated  using  piecewise  linear  finite  elements  and 
adaptive  h-refinement  on  an  initial  15  x  30  mesh.  The  color  red  denotes  a  high  concentra¬ 
tion  of  precursor  and  blue  denotes  a  low  concentration.  The  flow  enters  the  bottom  of  the 
reactor  and  exits  at  the  top.  The  fiber  is  heated  in  the  central  portion  of  the  reactor. 

For  low  values  of  Fr,  the  flow  circulates  and  less  precursor  exits  the  reactor  at  the 
top.  Larger  values  of  Fr  have  higher  speeds  and  less  precursor  gets  used  for  coating. 

The  Manufacturing  Group  of  Wright  Laboratories  has  been  given  a  copy  of  our 
software  and  has  been  using  it  to  conduct  parameter  studies  in  preparation  to  performing 
experiments.  They  will  also  use  it  to  design  a  neural  network  that  will  be  used  to  provide 
an  active  control  of  parameters  dming  fiber  coating. 

3.5.  Oxidation  of  Ceramic  Composites 

We  have  developed  a  model  to  study  the  oxidation  of  SiCISiC  composites  in  cracks. 
The  model  includes  solid-  and  gaseous  phase-reactions  and  viscous  momentum  equations; 
reactions  at  crack  surfaces;  and  bulk  and  Knudsen  diffusion.  Reactions  between  seven 
chemical  species  (silicon  carbide  5/C,  silicon  dioxide  5/0 2,  water  vapor  H2O,  oxygen  O2, 
carbon  monoxide  CO,  carbon  dioxide  CO 2,  and  hydrogen  arc  tracked.  The  model 
contains  features  similar  to  the  RVI  and  fiber  coating  models  of  the  two  previous  sections 
but  is  more  complex.  Nevertheless,  this  similarity  and  the  structured  design  of  the 
software  enabled  us  to  implement  the  without  excessive  coding. 

This  study  has  not  been  completed  but  work  is  continuing  and  results  will  be 
reported  [18].  A  preliminary  computation  appears  in  Figure  3.  The  figure  shows  half  of  a 
symmetric  idealized  crack  at  the  bottom  penetrating  a  5/C  matrix.  Surface  reactions 
liberate  O2  which  diffuses  into  the  matrix  and  oxidizes  the  5/C  (shown  in  red)  to  form  a 
5/0  2  layer  (shown  in  blue).  The  narrow  yellow  band  between  the  5/C  and  5/0  2  identifies 
the  reaction  zone.  CO  and  CO  2,  created  by  the  oxidizing  reaction,  diffuse  through  the 
solid  surface  and  escape  through  the  crack.  These  and  other  gaseous  effects  are  not 
shown  in  Figiu«  2. 

The  solid  is  extremely  viscous  at  the  operating  temperature  (1200  °C)  and  motion  of 
the  solid  is  not  apparent  at  the  recorded  times.  The  adaptive  mesh  is  concentrated  near 
the  reaction  zone  and  follows  the  front  as  it  progresses  through  the  matrix. 
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6.  Interactions  and  Transitions 

6.1.  Presentations 

1.  Flaherty  presented  an  invited  lecture  on  “Parallel  Adaptive  Finite  Element  Computa¬ 
tion”  at  the  Conf.  on  Application  of  High  Performance  Computing  in  Bioengineering, 
Pittsburg  Supercomputing  Center,  October,  21-22,  1994,  Pittsburg. 

2.  Haherty  presented  an  invited  lecture  “Adaptive  High-Order  Computational  Tech¬ 
niques  for  Singularly-Perturbed  Elliptic  and  Parabolic  Systems”  at  the  meeting  on 
Numerical  Methods  for  Singular  Perturbations,  Oberwolfach,  January  22-28,  1995. 

3.  Flaherty  presented  an  invited  lecture  “Adaptive  Computational  Methods  for  Compo¬ 
site  Materials  Processing”  Minisymp.  on  ODE  Methods  in  PDEs,  SCICADE  95,  Int 
Conf.  on  Sci.  Comput.  and  Diff.  ^ns.,  Stanford  University,  Palo  Alto,  March  28- 
April  1,  1995. 

4.  Flaherty  presented  an  invited  lecture  “Adaptive  Computational  Methods  for  Compo¬ 
site  Materials  Processing”  at  the  workshop  Applied  Mathematics:  Methods  and 
Applications,  Rensselaer  Polytechnic  Institute,  Troy,  March  24-25,  1995. 

5.  Flaherty  presented  a  seminar  “Parallel  Adaptive  Load-Balancing  Schemes  for  Three- 
Dimensional  Conservation  Laws”  at  New  York  University,  April  7,  1995,  New  York. 

6.  Adjerid  presented  a  lecture  “Modeling  and  Adaptive  Numerical  Solution  of  Fiber 
Coating  by  Chemical  Vapor  Deposition”  at  Aeromat  95,  ASM  Sixth  Int.  Aero.  Mats, 
and  Processes  Conf.,  Anaheim,  May  6-11,  1995. 

7.  Flaherty  presented  an  invited  lecture  “Adaptive  Method  of  Lines  Techniques  for  Par¬ 
abolic  Systems  with  Applications  to  Materials  Processing”  Workshop  on  the  Method 
of  Lines  for  Time  Dependent  Problems,  University  of  Kentucky,  Lexington,  May  31  - 
June  3,  1995. 

8.  Adjerid  presented  a  lecture  “A  Posteriori  Error  Estimation  for  Parabolic  Systems” 
Workshop  on  the  Method  of  Lines  for  Time  Dependent  Problems,  University  of  Ken¬ 
tucky,  Lexington,  May  31  -  June  3,  1995. 

9.  Haherty  presented  an  invited  lecture  “Adaptive  Computational  Methods  for  Compo¬ 
site  Materials  Processing”  at  the  symp.  on  Numerical  Methods  for  Industrial 
Manufacturing  Processes  1:  Composite  Materials,  Third  U.S.  Nat.  Congr.  of  Com¬ 
put.  Mech.,  Dallas,  June  12-14,  1995. 

10.  Adjerid  presented  an  invited  lecture  “Modeling  and  Adaptive  Numerical  Solution  of 
Fiber  Coating  by  Chemical  Vapor  Deposition”  at  the  symp.  on  Numerical  Methods 
for  Industrial  Manifacturing  Processes  1:  Composite  Materials,  Third  U.S.  Nat. 
Congr.  of  Comput.  Mech.,  Dallas,  June  12-14,  1995. 

11.  Haherty  presented  a  lecture  on  “Adaptive  and  Parallel  Computational  Techniques  in 
Materi^s  Science,”  at  the  1995  AFOSR  GranteesfContractors  Meeting  in  Computa¬ 
tional  and  Physical  Mathematics,  Albuquerque,  June  28-30,  1995. 

12.  Haherty  presented  an  invited  lecture  “Adaptive  Methods  for  Parabolic  Partial 
Differential  Equations  with  Applications  to  Shear  Band  Formation”  at  ICES  95,  Int 
Conf.  on  Comput  Engng.  Sci.,  Mauna  Lani,  July  30  -  Aug.  3,  1995. 
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13.  Flaherty  presented  an  invited  lecture  “Adaptive  and  Parallel  Computational  Tech¬ 
niques,”  at  the  IBM  Symposium  on  Parallel  Computation,  Supercomputing  ’95,  San 
Diego,  December  4-8,  1995. 

14.  Flaherty  presented  a  seminar  on  “Adaptive  and  Parallel  Computational  Techniques 
with  Applications  in  Materials  Science,”  Air  Products,  Inc.,  Allentown,  January  23, 
1996. 

15.  Flaherty  presented  the  Schiesser  Lecture  in  Interdisciplinary  Science  and  Engineering 
on  “Adaptive  Method  of  Lines  Techniques  for  Parabolic  Systems  with  Applications 
to  Materials  Processing,”  Lehigh  University,  January,  24,  1996. 

16.  Flaherty  presented  a  seminar  on  “Adaptive  and  Parallel  Computational  Techniques 
for  Partial  Differential  Equations,”  U.S.  Army  High  Performance  Computation 
Center,  University  of  Minnesota,  February  8  1996. 

17.  Flaherty  presented  a  seminar  on  “Adaptive  and  Parallel  Computational  Techniques 
for  Partial  Differential  Equations,”  Theory  Center,  Cornell  University,  March  4, 
1996. 

18.  Flaherty  presented  a  seminar  on  “Adaptive  and  Parallel  Computational  Techniques 
for  Parti^  Differential  Equations,”  Department  of  Applied  Mathematics,  Brown 
University,  March  8,  1996. 

19.  Flaherty  gave  a  talk  on  “Adaptive  and  Parallel  Computation  for  Materials  Process¬ 
ing”  at  the  AFOSR  Grantees  and  Contractors’  Meeting:  Computational  and  Physical 
Mathematics,  Wright  Laboratories,  Wright  Patterson  Air  Force  Base,  June  24-26, 
1996. 

20.  Flaherty  gave  a  keynote  lecture  on  “High-Order  Adaptive  Computational  Methods 
With  Applications  in  Materials  Science”  at  the  Conference  on  Grid  Adaptation  in 
Computational  PDEs:  Theory  and  Applications,  University  of  Edinburgh,  Edinburgh, 
July  1-5,  1996. 

6.2.  Advisory  Functions  to  Laboratories  and  Agencies 

1.  Adjerid,  Flaherty,  and  Hudson,  visited  the  Manufacturing  group  at  Wright  Labora¬ 
tories,  Wright-Patterson  Air  Force  Base,  on  December  12,  1994.  We  met  with 
Steven  Le  Clair,  James  Malas,  John  Busbee,  et  al.  and  discussed  om  research  on  the 
coating  of  sapphire  fibers  in  relation  to  their  experimental  work  on  the  same  topic. 
Bruce  Webster  of  Centric  Engineering  was  also  in  attendance. 

2.  Adjerid  and  Flaherty  collaborated  with  Robert  Ferencz,  Thomas  Hughes,  and  Bruce 
Webster  of  Centric  Engineering  on  modeling  and  the  numerical  solution  of  fiber  coat¬ 
ing  processes.  Ferencz  visited  Rensselaer  on  July  12-14,  1995  and  Flaherty  visited 
Centric  on  March  30  and  July  28,  1995. 

3.  Adjerid  and  Flaherty  met  with  John  Gamier  of  Du  Pont  Lanxide  to  discuss  our 
mutual  interests  in  processes  to  coat  ceramic  fibers,  June  1,  1995.  Lanxide  is  provid¬ 
ing  silicon  carbide/silicon  carbide  composite  samples  to  Rensselaer  researchers  for 
mechanical  testing.  These  composites  contain  coated  fiber  tows  and  their  data  could 
provide  an  excellent  platform  for  testing  our  models  and  software. 
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4.  Richard  K.  Everett  of  the  Naval  Research  Laboratory  has  expressed  an  interest  in  col¬ 
laborating  with  us  on  an  investigation  involving  the  reaction  and  diffusion  between 
nickel  and  aluminum  films  in  multilayers.  Adjerid  and  Raherty  met  with  Everett  on 
June  1,  1995  for  an  initial  discussion  of  this  possibility. 

5.  Flaherty  visited  Karen  Devine,  Andrew  Salinger,  and  John  Shadid  of  Sandia  National 
Laboratories,  June  28,  1995  to  discuss  our  parallel  adaptive  software  for  solving 
CVD  problems. 

6.  Adjerid  and  Flaherty  visited  John  Jones  of  Wright  Laboratories  on  June  24,  1996  and 
discussed  our  research  on  fiber  coating. 

6.3.  Transitions 

1.  Adjerid  and  Flaherty  (in  collaboration  with  Hudson  and  Shephard  of  Rensselaer  and 
Robert  Ferencz,  Thomas  Hughes,  Bruce  Webster  of  Centric  Engineering)  developed  a 
mathematical  model  and  computer  software  for  analyzing  axisymmetric  fiber  coating 
processes  by  CVD.  John  Jones  of  Wright  Laboratories  has  a  copy  of  this  software 
and  has  been  using  it  for  parameter  studies  and  to  create  a  neural  network  model  of 
the  parameter  space  that  will  be  used  to  control  their  experiments  on  the  coating  of 
sapphire  fibers. 

2.  Our  fiber-coating  model  and  adaptive  software  will  be  encorporated  into  Centric 
Engineering’s  SPECTRUM  finite  element  software  for  solving  stress  and  fluid-flow 
problems. 


